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Abstract 

We present a theory of neutrino interactions with nuclei aimed at the descrip- 
tion of the partial cross-sections, namely quasi-elastic and multi-nucleon emission, 
coherent and incoherent single pion production. For this purpose, we use the the- 
ory of nuclear responses treated in the random phase approximation, which allows 
a unified description of these channels. It is particularly suited for the coherent 
pion production where collective effects are important whereas they are moderate 
in the other channels. We also study the evolution of the neutrino cross-sections 
with the mass number from carbon to calcium. We compare our approach to the 
available neutrino experimental data on carbon. We put a particular emphasis on 
the multi-nucleon channel, which at present is not easily distinguishable from the 
quasi-elastic events. This component turns out to be quite relevant for the interpre- 
tation of experiments (K2K, MiniBooNE, SciBooNE). It can account in particular 
for the unexpected behavior of the quasi-elastic cross-section. 
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1 Introduction 

Neutrino physics has undergone a spectacular development in the last decade, follow- 
ing the discovery of neutrino oscillations first revealed by the anomaly of atmospheric 
neutrinos [l|. A number of results on the interaction of neutrinos with matter are now 
available. Neutrino detectors do not usually consist of pure hydrogen but they involve 
complex nuclei for instance 12 C, as in SciBar [2i], where the molecule CsH 8 is involved, 
or in MiniBooNE 0| which uses the mineral oil CH 2 . Heavier nuclei are also under 
consideration for instance in the liquid argon chamber planned for T2K 0, 0|- A num- 
ber of results have been obtained, for neutral or charged current (K2K, MiniBooNE, 
SciBooNE) on quasi elastic processes or coherent and incoherent single pion production 
@, 0, E3, El, H Q Q H 0, 03- The first question is then if our present un- 
derstandin g of neutrino interactions with matter can re pro duce the available data. Many 



works [JJ, bj, |20j, |2l|, |22|, |2J, M, \M M, |27j, 12$ |29|, |30|, l3l|, |3J, 1MJ34JJ35 



been devoted to this problem, using various theoretical approaches [37|, |38|, |39j, |40j, 141 



36| have 
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42|, |43|, m , m m m m, m Lso, y, m, m m, m \m m m m, m M, lea m m 



In this article we will explore such interactions using the theory of the nuclear response 
treated in the random phase approximation (RPA) in the quasi-elastic and Delta reso- 
nance region including also two and three nucleon knock-out. The formalism is the same 
as the one used by Marteau 
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in his work on the v— 16 interaction. The merit of 
this approach is that, although perfectible in several ways, it describes in a unique frame 
several final state channels. This technique has been successful in a number of problems 
involving either weakly interacting probes such as (e,e') scattering or strongly interacting 
ones such as pion scattering or ( 3 He,T) charge exchange reaction [65]. We give the cross- 
sections for pion production, coherent or incoherent, and nucleon knock-out, for neutral 
or charged currents. We restrict to single pion production ignoring two-pion production 
processes which, for real photons, lead to a sizable part of the photo-absorption cross- 
section at energies larger than the Delta resonance, above ~ 500 MeV. Our treatment 
should thus underestimate the cross-section when multi-pion production starts to show 
up. Our work ignores as well the meson exchange effects which play a non negligible role 



66j, |67fl. We only take into account the exchange effect in the time component of the axial 
current, which is known to be important (68|. For single pion production we assume that 
the dominant production mechanism is via the Delta resonance, ignoring the other reso- 
nance excitations, which also limits the energy for the validity of our approach. Beyond 
quasi-elastic processes and single pion production via Delta excitation we also incorporate 
several nucleon knock-out through two-particle-two-hole (2p — 2h) and 3p — 3h excitations. 
These will play a crucial role in the comparison with data involving quasi-elastic events. 

Among the aims of this work there is the exploration of the evolution of the neutrino- 
nucleus interaction as the mass number of the nucleus goes from the carbon region to the 
region of 40 Ca. This investigation is motivated by the project of a liquid argon chamber in 
the T2K experiment which raises the question if one keeps control of the understanding 
of the interaction of neutrinos with matter by going to a medium-weight nucleus such 
as 40 Ar. In order to single out the evolutions linked to the nuclear size we have chosen 
as element of comparison an isoscalar nucleus in the 40 Ar region, namely 40 Ca. For the 
coherent process which per nucleon fades away in heavy nuclei, the evolution is relatively 
rapid but should remain under control as our theory is particularly well adapted to this 
channel. The other exclusive channels, in particular the incoherent pion production, are 
sensitive to final state interaction not automatically included in our approach. This leaves 
some uncertainty in the evolution between the mass 12 and the mass 40 region for this 
channel. 

Our article is organized as follows: Section 2 introduces the formalism of the response 
functions treated in the random phase approximation (RPA). Section 3 discusses the 
various final state channels. In Sec. 4 we compare these predictions with the available 
data. In Sec. 5 we provide a summary and conclusion of the present work. 



2 



2 Formalism 



The double differential cross-section for the reaction vi (yj) + A — > I + X is given 
by 

d 2 a _G 2 F cos 2 9 c k' 2 



dVL k ,dk' 327T 2 k k' 

where Gp is the weak coupling constant, 9 C the Cabbibo angle, k and k' the initial and final 
lepton momenta, and T the invariant amplitude given by the contraction of the leptonic 
L and hadronic W tensors. Their expressions are given in Appendix [A] In order to 
illustrate how the various response functions enter and to introduce the variables, we give 
below a simplified expression, which in particular ignores the lepton mass contribution 
and assumes zero A width. We stress however that in the actual calculations the full 
formulas of Appendix [A] which do not make these simplifications, have been used. The 
simplified double differential cross-section reads 
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(2) 



where % = k^ — k' ll = (u, q) is the four momentum transferred to the nucleus, 6 the scat- 
tering angle, M& (M n ) the Delta (nucleon) mass. The electric, magnetic and axial form 
factors are taken in the standard dipole parameterization with the following normaliza- 
tions: Gp(0) = 1.0, Gm(0) = 4.71 and Ga(0) = 1.255. The corresponding cut-off param- 
eters are My = 0.84 GeV/c 2 for the electric and magnetic terms and Ma = 1.032 GeV/c 2 
for the axial one. The plus (minus) sign in Eq. ([2|) stands for the neutrino (antineutrino) 
case. A similar expression applies to the process : v\ {pi)+A — > ui (z/;)+X, which involves 
neutral currents. The various responses R appearing in Eq.((2j) are defined according to 



A A 

K P ' =E<"lE ^l°)H E 6 ^ ~ E » + E 0)- (3) 

n j=l k=l 

The upper indices (P, P') refer to the type of particle (N or A) at the vertices that couples 
to the external probe. The corresponding operators have the following forms: 

0»(j) = rf, (aj.qjrf, (<r,xg)Vf, 

for a = r, err(L), ar(T), and 
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for a = crr(L), <jt(T). We have thus defined the inclusive isospin (R T ), spin-isospin 
longitudinal (R aT ^) and spin-isospin transverse (R aT (T)) nuclear response functions (the 
longitudinal and transverse character of these last two responses refers to the direction 
of the spin operator with respect to the direction of the transferred momentum). The 
operators S and T are the usual 1/2 to 3/2 transition operators in the spin and isospin 
space. We have assumed the existence of a scaling law between the nucleon and Delta 



magnetic and axial form factors 69 



G*m/Gm — G* a /Ga — /*//, 

where /* (/) is the n N A (ir N N) coupling constant. For a matter of convenience, we 
have incorporated the scaling factor /*// = 2.2 into the responses. 

The presence of the spin-isospin longitudinal coupling is a distinct feature of neu- 
trino interaction as compared to inelastic electron scattering. For instance coherent pion 
production, present in v interactions is partly suppressed in (e, e') scattering due to the 
purely transverse spin coupling of the exchanged photon. Inclusive electron scattering is 
nevertheless useful as a test for the transverse response 64|. The response functions are 



related to the imaginary part of the corresponding full polarization propagators 

R(u,q) = --lm[U(u,q,q)}, (4) 

71 

where V is the nuclear volume such that Vp = A. They are calculated within a RPA 
(random phase) ring approximation starting from "bare" propagators (meaning that the 
nuclear correlations are switched off) . The word bare here does not imply that the corre- 
sponding response is free of many-body effects, as described in the following. The "bare" 
polarization propagator is illustrated by some of its components in Fig. Q] where the wig- 
gled lines represent the external probe, the full lines correspond to the propagation of a 
nucleon (or a hole), the double lines to the propagation of a Delta and the dashed lines 
to an effective interaction between nucleons and/or Deltas. 

The dotted lines in Fig. CD indicate, in each of the channels introduced previously 
(NN,NA or AA), which intermediate state is placed on-shell. It follows that the bare 
response is the sum of the following partial components: NN : quasi-elastic (as described 
by the standard Lindhard function), NN : 2p - 2h, NA and AiV : 2p - 2h, AA : ttN, 
AA : 2p- 2h, AA : 3p - 3h. 

Notice that the graphs shown in Fig. [T] do not exhaust all the possibilities for the bare 
propagator. For instance the distortion of the pion emitted by the A is not explicitly 
shown, although it will be included in our evaluation through the modification of the A 
width in the nuclear medium. But the type of final states that we consider is limited to the 
previous list. Thus, in the bare case, through the introduction of the partial polarization 
propagators illustrated by the Feynman graphs of Fig. CO the inclusive expression of Eq. 
j2j) provides an access to the exclusive ones, with specific final states. 

For the actual evaluation of the bare response, i.e., the imaginary piece of the bare 
propagator, some of the graphs of Fig. Q] amount to a modification of the Delta width 
in the medium. We take into account this modification through the parameterization of 
the in-medium Delta width of Oset and Salcedo which leads to a good description 
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Figure 1: Feynman graphs of the partial polarization propagators: NN quasi-elastic (1), 
NN (2p-2h) (2), NA (2p-2h) (3), AN (2p-2h) (3'), AA (irN) (4), AA (2p-2h) (5), AA 
(3p-3h) (6). The wiggled lines represent the external probe, the full lines correspond to 
the propagation of a nucleon (or a hole), the double lines to the propagation of a Delta 
and the dashed lines to an effective interaction between nucleons and/or Deltas. The 
dotted lines show which particles are placed on-shell. 



of pion-nuclear reactions. The authors split the A width into different decay channel 
contributions : the A — > n N, which is modified by the Pauli blocking of the nucleon 
and the distortion of the pion. Moreover, in the nuclear medium, new decay channels are 
possible: the two-body (2p — 2h) and three-body (3p — 3h) absorption channels which 
they also incorporate. They give a parametrization for the inclusion of these effects, 
both in the case of pion interaction with nuclei and for the photo-production process. 
We have used their parametrization in spite of the fact that in neutrino interaction the 
intermediate boson has a space-like character. An explicit evaluation of the corresponding 
contributions in the kinematical situation of neutrino scattering is desirable. There exist 
also other 2p — 2h contributions which are not reducible to a modification of the Delta 
width. We include them, as in the work of Marteau (38|, following the method of Delorme 



and Guichon [71fl who perform an extrapolation of the calculations of Ref. [72j on the 
2p — 2h absorption of pions at threshold. For the last contribution only the imaginary 
part of the corresponding propagator is incorporated. The explicit expressions are given 
in Appendix lB.il It turns out that for neutrino interaction it is the dominant contribution 
to the 2p — 2h final state channel, as will be illustrated later. This piece of the cross- 
section is subject to some uncertainty as this parametrization has not been constrained 
by specific experimental tests. This point will be discussed in more detail in Sees. 13.31 
and ELS 

The "bare" polarization propagator is density dependent. In a finite system, II (u, q, q'), 
it is non-diagonal in momentum space. In order to account for the finite size effects we 
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Figure 2: Bare response for 12 C at q = 300 MeV/c as a function of the energy transfer 
with its different components, quasi-elastic and pion emission (left panel), 2p — 2h and 
3p — 3h (right panel). 



evaluate it in a semi-classical approximation where it can be cast in the form 

n V, q, q') = J dr e<*-*> r n° (u, l -{ q + q >) , r) . (5) 

In practice we use a local density approximation, 

tf( w ,i±*r)=nS, w («,2±*), (6, 

where the local Fermi momentum kp(r) is related to the experimental nuclear density 
through : k F (r) = (3/2 n 2 p(r)) 1 / 3 . The density profiles of the various nuclei considered 
are taken from the Sum-of-Gaussians nuclear charge density distribution parameters ac- 
cording to Ref. The corresponding bare response for 12 C at q = 300 MeV/c as a 
function of the energy transfer is illustrated in Fig J2] with its different components, quasi- 
elastic, pion emission, 2p — 2h and 3p — 3h. In all figures the responses incorporate the 
multiplicative spin-isospin factor. 

Turning to the Random Phase Approximation, as the semi-classical approximation is 
not suited to evaluate the collective effects, we have used the previous bare polarization 
propagator IT as an input in a full quantum mechanical resolution of the RPA equations 
in the ring approximation. The introduction of the RPA correlations amounts to solving 
integral equations which have the generic form: 

n = n° + n° vu, (7) 

where V denotes the effective interaction between particle-hole excitations. Its diagram- 
matic representation is given in FigEl Some detailed expressions are given in Appendix 
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IB.2L In the spin-isospin channel the RPA equations couple the L and T or the iV and 
A components of the polarization propagators. The effective interaction relevant in the 
isospin and spin-isospin channels is the crucial ingredient for determining the importance 
of the RPA effects. We use the parametrization in terms of it, p and contact pieces: 



V NN = (f + K- + V p + V gl ) n.r 2 

V NA = (K + V p + V g .) t x .T\ 

V an = (V w + V P + V g ,) Ti.t 2 

Vaa = (K + V P + V g> ) T X .T\. (8) 



For instance, in the A^A" case one has : 



> 2 o 
9r \ „2 <1 



\2M N J lu 2 - q 2 - m 2 

(\ 2 2 
■Sr) c p f p -2—^2 i<nxq<72xq 
2M N J ' p to 2 — q 2 — m 2 



where g' is the Landau-Migdal parameter and C p = 1.5. Here F n (q) = (A^ — m^)/(A^ — q 2 ) 
and F p (q) = (A 2 p — m 2 p )/(A 2 p — q 2 ) are the pion-nucleon and rho-nucleon form factors, 
with A„. = 1 GeV and A p = 1.5 GeV. For the Landau-Migdal parameter /', we take 
/' = 0.6. As for the spin-isospin parameters g' we use the information of the spin-isospin 



phenomenology [7j|, with a consensus for a larger value of g' NN = 0.7; for the other 
parameters we take g' NA = g' AA = 0.5. 

The separation between the specific channels is less straightforward in the RPA case 
than in the bare one. Indications can be obtained with the following method, introduced 



in Ref. [38||. The imaginary part of II can be written (again generically) as : 

Imn = |n| 2 ImV + |1 + UV\ 2 mill . (10) 

It separates into two terms. The first term on the r.h.s. of Eq. (PT0l) . IllpImV^, is absent 
when the effective interaction is switched off. In the domain of energy considered it is the 
imaginary part of the pion exchange potential V w which plays the major role. This process 
thus represents the coherent pion production, i.e., the emission of an on-shell pion, the 
nucleus remaining in its ground state. This is illustrated in Fig. [H in which the hatched 
rings represents the RPA polarization propagator. The second term on the r.h.s. of Eq. 
(fTO]) . proportional to the bare polarization propagator Imll , reflects the type of final 

state already mentioned for the imaginary part of n : NN,tcN, The factor in front, 

|1 + II V| 2 , embodies the modification of the exclusive bare responses by the collective 
effects. We point out however that final state interactions are not incorporated in this 
description. For instance a pion produced in the decay of the Delta resonance can be 
absorbed on its way out leading to a multi-nucleon emission process. Thus the second 
term in Eq. (ITOll is adequate for the sum of the incoherent pion production and the 
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Figure 3: Diagrammatic representation of the RPA polarization propagator. The white 
bubble is the free p-h propagator while the black is the full RPA one. 




Figure 4: Diagrammatic representation of the coherent process. The dotted line indicates 
that pion is placed on-shell. 



multi-nucleon knock-out channels, but not for each channel individually. The separation 
between these two channels from the type of final state is approximate for light nuclei such 
as 12 C. In heavier nuclei it overestimates the incoherent pion channel, underestimating 
the multi-nucleon one. We will illustrate this fact in the scattering of physical pions. 

Having established the formalism, we are now ready to evaluate the cross-sections 
in the various partial channels. In the actual numerical calculation we have limited the 
energy transfer to uj — 1 GeV as our approach becomes insufficient for a larger energy 



transfer. The center-of-mass correction for the n-N system gcM = i+^m UM i s made by 
dividing the bare responses by a factor r 2 = (1 + u/M) 2 . The components of the neutrino 
cross-section that does not involve the momentum q at the two ends of the RPA chain 
are obtained by an overall multiplication by the factor r 2 . Interference terms with one 
momentum are multiplied by r. 
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Figure 5: Longitudinal and transverse total responses of 12 C at fixed q = 300 MeV/c as a 
function of u. The coherent part of the responses is also shown. 

3 Results 

3.1 Coherent cross-section 

Several types of responses enter the total neutrino cross-section, isovector, spin-isospin: 
transverse or longitudinal. The last quantity is naturally associated with the coherent 
process, since it has the same coupling as the pion. The production by a transverse spin 
coupling requires a transverse-longitudinal conversion which is partly suppressed. This 
difference is illustrated in Fig. [5] where the total responses, longitudinal and transverse, 
of 12 C are displayed as a function of the energy transferred to the nuclear system for 
a fixed three-momentum q= 300 MeV/c. The coherent component, much larger in the 
longitudinal case, is also shown. 

Figure [6] illustrates the evolution with the nuclear size of the coherent part of the 
longitudinal response per nucleon as a function of the energy at fixed momentum for 
some nuclei, 12 C, 16 0, 40 Ca and also for a fictitious piece of isospin symmetric nuclear 
matter with the density profile of lead. Two features emerge, the first one is that its 
magnitude decreases in "lead", as expected : the coherent response per nucleon vanishes 
in nuclear matter when the polarization propagators become diagonal in momentum space. 
The second is that the coherent response is not peaked at the energy ui n = (q 2 + m 2 ) 1 / 2 
where the mismatch between the incident energy and that of the physical outgoing pion 
is smallest. Instead it is reshaped by the collective features of the longitudinal response 
with the appearance of two collective branches on each side of the pion line. This is more 
apparent in the case of the (fictitious) lead. 

As a test of our description of the coherent responses we have investigated the elastic 
scattering of pions on nuclei in the Delta region, related to the coherent part of the 
spin-isospin longitudinal response through : 
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Figure 6: Evolution with the mass number of the coherent longitudinal response per 
nucleon at fixed q = 300 MeV/c as a function of u>. The arrow indicates the energy for 
on-shell pion. 
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Figure 7: 7r- 12 C elastic cross-section as a function of pion energy. 
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Figure 8: Charged and neutral current coherent pion production differential cross-section 
off 12 C versus pion kinetic energy for several energies. 
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Figure 9: Total CC (divided by 2) and NC z/^-induced coherent pion production cross- 
sections in 12 C function of neutrino energy. 
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= Uj*A 2 R c L oh (u, o (ii) 

where q\ = u 2 — ml and Rf h refers to the coherent part of the longitudinal response. The 
resulting cross-section in the case of 12 C is shown in Figj7| together with the experimental 
points from Ref.[76]. The agreement with data is satisfactory. A similar accuracy can be 
expected for the coherent response which enters the neutrino cross-section, at least in the 
energy region for the produced pion where we have tested our model (i.e., between u ~ 220 
and ~ 450MeV r ). The elastic cross-section which depends on the longitudinal response 
is particularly sensitive to collective effects in this channel known to be important. The 
replacement of the bare response by the RPA one leads to a different energy behavior, 
the collective effects producing a softening of the response, characteristic of the collective 
nature of the longitudinal channel. 

Figured] displays our evaluations of the neutrino coherent cross-section on 12 C as a function 
of the pion kinetic energy, both for charged and neutral current, for several neutrino 
incident energies. The resulting total coherent cross-sections are displayed in Figj9] 
The suppression of the meson exchange correction in the time component of the axial 
current, G* A — > Ga, produces a moderate ~ 10% increase of the cross-section. 
The data available on the coherent production by neutrino concern its ratio to the total 
cross-section and to the total pion production. We will then postpone the comparison 
with experimental data after the discussion of the various other channels. 



3.1.1 Adler's theorem 

In the forward direction where q = u and for vanishing lepton mass, only the spin longitu- 
dinal response contribution survives. As it also enters in pion scattering, it is possible to 
relate the forward neutrino cross-section to the cross-section of physical pions, apart from 
a difference in kinematics: q = to (soft pions) for neutrinos, instead of q = q w = ■sJuS 1 — m\ 
for physical pions. This difference becomes less relevant at large energies. This is the con- 



tent of Adler's theorem [73]. The coherent channel, which is completely dominated by 
the longitudinal response, offers the best application of this theorem, while for the other 
channels the transverse component, which bears no relation to pion scattering, quickly 
takes over as soon as one moves away from the forward direction. This theorem has been 



used in the approach of Refs. [18| [29|] [3J| to evaluate the coherent neutrino-nucleus 
cross-section. This is not our aim here. We want to illustrate the link between the for- 
ward direction coherent neutrino cross-section and the elastic pion-nucleus one. For the 
coherent cross-section Adler's relation writes 

d 2 o \ coh G 2 F cos 2 9 C f2 E u -uj las 



dil duj ) g=Q 7r 3 77 UJ 



tt— (12) 



where f n = 93.2 MeV is the neutral pion decay constant. Introducing the experimental 
values for the elastic cross-section taken from Ref.]76] we obtain the points shown in Fig. 
[TOl together with our predicted curve. The agreement is rather good. It deteriorates at 
small energies when the kinematical difference between soft and physical pions becomes 
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Figure 10: v e - 12 C coherent cross-section in the forward direction. Continuous line: our 
result. Circles: deduced, according to Adler's relation of Eq. (fT2l) . from the experimental 
values for the elastic cross-section taken from Ref. [76| . Stars: introducing into the r.h.s. 
of Eq. (fT2j) the multiplicative factor — . 



substantial. A natural correction can be performed with the introduction into the r.h.s. 
of Eq. (fT2l) of a multiplicative factor as suggested by the relation of Eq.JTT]) between 
Rl and a elas \ The corresponding corrected points are also shown in Fig. [TOl extending 
somewhat the region of agreement. The use of the Adler relation becomes problematic at 
energies near threshold. For small neutrino energy (E u < 0.5 GeV) this region has more 
weight in the total coherent cross-section. 

The Adler relation thus provides a good test for our evaluation of coherent neutrino 
cross-section in the forward direction. We believe that the extrapolation to the non- 
forward direction as performed in our model should be under control. 

3.2 Pion-nucleus cross-sections 

The various partial cross-sections for physical pions on nuclei constitute a precious piece of 
information. Elastic cross-section has already been introduced as a test for the coherent 
cross-section. The total cross-section for pions on the nuclei is given by an expression 
similar to Eq.lfTTT) with the full polarization propagator replacing the coherent piece 



The corresponding cross-section is displayed in FigfTT] together with the experimental 
points. We will show that in the same way the inelastic cross-section provides some 




(13) 
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information on the incoherent pion production by neutrinos and the absorptive cross- 
section on the multi-nucleon channels. Figure [TT] displays the various partial channels 
(but the elastic one, previously shown) which contribute to the 7r + cross-section on 12 C, 
namely the inelastic pion scattering channel (which is the incoherent scattering with a tt + 
in the final state) and the absorptive one. We also display the sum of the incoherent pion 
(including charge exchange) and true absorption (multi-nucleon channels) cross-sections. 
The experimental points are taken from Ashery et al. (76| . To reduce the clutter, we 
have not explicitly plotted the charge-exchange cross-section which, in our approach, is 
one fifth of the inelastic tt + cross-section and is consistent with the experimental data. 
While the elastic cross-section was well reproduced, our approach overestimates the 7r + 
inelastic channel in the peak region and largely underestimates the absorptive channel. 
We attribute this deficiency to the absence of pion final state interaction as the pion can 
be reabsorbed on its way out the nucleus. It can also undergo charge exchange process 
but this is a smaller effect. As a counterpart the absorptive multi-nucleon production 
is underestimated, as is apparent in FigfTU The sum of the two channels is instead 
reasonably well reproduced in the peak region. 

These limitations also affect the incoherent neutrino-nucleus cross-section but we stress 
that, in contradistinction, our description for the coherent channel automatically contains 
the final state interactions and no further correction is needed. The total neutrino cross- 
section is also obviously not affected. With the information on the pion energy spectrum 
in neutrino interactions (that our calculation does not provide) it would be possible to 
estimate at each energy an attenuation factor for the incoherent neutrino production 
from the difference between our calculation and inelastic data for physical pions. For 
instance, for 12 C at E v — 1 GeV, a rough evaluation of the overall correction for the 
incoherent production cross-section with the information on the pion spectrum [l8| results 
in a moderate reduction of ~ 15%. A similar attenuation was found in oxygen at E v = 500 



MeV and E v = 750 MeV [43|. A larger correction is obviously expected for calcium. 



3.3 Quasi-elastic and multi-nucleon channels 

The quasi-elastic (QE) channel corresponds to a single nucleon knock-out. In the quasi- 
elastic process the space-like character is pronounced as the quasi-elastic peak occurs at 
uj ~ q 2 /(2M7v), hence the distribution in Q 2 = q 2 — to 2 is rather broad @. At zero order 
only R NN contributes to this channel. In the RPA chain instead R NA and R AA also 
participate. For instance the lowest order contribution of R NA is illustrated in Fig. Q21 

In contrast to the coherent channel, the quasi-elastic one is totally dominated by the 
transverse response. The longitudinal contribution is suppressed by a cancellation between 



the space and time components of the axial current, as observed by Marteau 38] and 
shown in Appendix IA.ll for vanishing lepton mass and neglecting the Fermi momentum. 
Numerically its contribution is indeed very small. We have tested our semi-classical 
approximation on the bare QE y P - 12 C cross-section through a comparison with the one 



obtained by Martini et al. [531 ] in the continuum shell model where the mean field is 



produced b y a Woods-Saxon well. Our result is very similar in shape and magnitude to 



the one of [53j but for a displacement in energy of 27 MeV. This reflects the inclusion 



of the nucleon separation energy in the continuum shell model, which is ignored in our 
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Figure 11: Partial and total ix - 12 C cross-sections. 




,ure 12: Lowest order contribution of R to the quasi-elastic channel. 
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approximation. 

The quasi-elastic cross-section is displayed in Fig. [13] as a function of the energy 
transfer for neutrino energy E v — 1 GeV, both in the bare case and in the RPA one. 
The RPA influence produces a reduction, as expected from the repulsive character of 
the particle-hole interaction, which prevails in the transverse channel. This reduction is 
mostly due to the interference term R NA which is negative (Lorentz-Lorenz effect (79|). 

The total quasi-elastic charged current and neutral current cross-section are plotted 
in FigsHH and [[5] as a function of the neutrino energy. In Figs. [EH El an d CES we also 
display the sum of the two- and three-nucleon knock-out cross-sections, which represents 
a sizable fraction of the quasi-elastic one. Singling out the genuine quasi-elastic process 
requires the insurance that no more than one proton is ejected. This question will appear 
in the comparison with data. Among the various contributions to the multi-nucleon 
channel the ones which do not reduce to a modification of the A width are dominant. The 
accumulation of 2p—2h strength at low energy is an artifact of the simplified extrapolation 
that we use in this channel. In Section l4~3~l this point is discussed in more detail and another 
method for the parametrization, with an explicit momentum dependence, is introduced. 
It modifies the u dependence of spreading the strength over a larger energy region, 
but does not substantially affect the energy integrated cross-section. 

Coming now to the evolution of these channels between 12 C and 40 Ca we compare the 
corresponding RPA differential cross-sections per neutron for the two nuclei in Fig. HH 
One can see that the evolution of this quantity with the mass number is quite weak in 
the QE case. It is also weak in the multi-nucleon channel although it should increase 
faster with density than the quasi-elastic one. However, between a light system such as 
12 C and 40 Ca the evolution is moderate. Only in the case of deuteron one expects the 
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Figure 14: CC - 12 C QE, multi-nucleon [np — nh), incoherent pion emission and total 
cross-section as a function of neutrino energy. 



multi-nucleon knock-out to be appreciably smaller in view of the loose binding of the 
system. 

3.4 Incoherent pion emission 

The pion arises from the pionic decay of the Delta leaving the nucleus in a p — h excited 
state. For the nuclei that we consider this cross-section is much larger than the coherent 
one. As compared to a free nucleon the emission probability is already appreciably reduced 
in the bare case by the change in the Delta width. Moreover the RPA effects, which are 
moderate, also tend to a small reduction. The reduction due to the modification of the 
Delta width has a counterpart in the presence of a component of multi-nucleon knock- 
out. Charged current and neutral current cross-sections for incoherent pion emission for 
all possibles charges are represented in FiglH] and [15] as a function of neutrino energy. 
Moreover these figures summarize all previous results for the other channels and also give 
the total cross-sections. 

On the other hand, FigJT6l compares the neutrino differential cross-section per neutron 
in the various channels as a function of the energy transfer, u, for the cases of 12 C and 
40 Ca and for a neutrino energy E v — 1 GeV. The two sets of curves are very similar. We 
can conclude that, at the level of our approximation, i.e., without final state interaction, 
it is possible to extrapolate smoothly from 12 C to the region of 40 Ar. Only the coherent 
cross-section presents a significant variation, illustrated in FigJT7l 
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cross-section as a function of neutrino energy. 
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Figure 16: RPA differential CC cross-sections per neutron in the different channels for 
12 C (same convention-line as Figs JT3l|T4lfT5|) . and 40 Ca (dotted lines). 
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Figure 17: Differential CC cross-sections per neutron in the coherent channel for 12 C 
(continuous line), and 40 Ca (dotted line). 



4 Comparison with data 
4.1 Coherent pion production 

Experimental data concern ratios between different cross-sections. The first indication of 
coherent pion production by neutral current was given by MiniBooNE |l0|, that found the 
ratio of coherent to total 7T° production to be 0.195 ± 0.011 ± 0.025. In this experiment 
the neutrino flux is spread in energy with a peak at ~ 700MeV [13]. Our approach 
leads to a lower number, namely 0.06, which is difficult to reconcile with experimental 



data, a problem that other groups also face. It has been suggested in Ref. 27] that 



MiniBooNE, which uses Rein-Sehgal model [80j for data analysis, possibly overestimates 



the 7T° coherent cross-section. In a preliminary report [8l| the experimental value given 
for this cross-section is 7.7±1.6±3.6 10~ 40 cm 2 . Our result for this cross-section averaged 
on the MiniBooNE flux (la ]. 2.8 10~ 40 cm 2 , is compatible with the experiment in view of 
the large experimental errors. 

On the other hand for charged current, two experimental groups have given upper 
limits for the ratio of coherent pion production to the total cross-section. The K2K 
collaboration gives a limit of 0.60 10~ 2 averaged over a neutrino flux with a mean energy 
of 1.3 GeV 0]. More recently, the SciBooNE collaboration found for the same quantity 



0.67 10 -2 at neutrino energy of 1.1 GeV [l2| and 1.36 10~ 2 at neutrino energy of 2.2 GeV. 
We report in FigH8] our prediction for this quantity. Since our approach is appropriate 
for a limited neutrino energy range we keep in the comparison only the lowest energy 
SciBooNE point. Our curve is just compatible with the experimental bound. 
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Figure 18: Ratio of the z/^-induced charged current coherent 7r + production to total cross- 
section as a function of neutrino energy. 



4.2 Total pion production 

Another measured quantity is the ratio of 7r + production to quasi-elastic cross-section 
for charged current. The MiniBooNE collaboration has used a CH 2 target. In order to 
compare with ANL [82] and K2K [IH data, they presented the results with an isoscalar 
rescaling correction [lj]. The issue of pion loss by final state interaction, which is not 
incorporated in our description, has also been taken into account by MiniBooNE who 
corrects data for this effect. We can thus compare our n + over quasi-elastic ratio (solid 
line in the upper panel of Fig fl9]) to the final-state-interaction-corrected MiniBooNE re- 
sults. Our curve incorporates the small coherent cross-section; the incoherent pion one 
is multiplied by the isospin factor 5/6 to single out n + contribution. Our curve is fully 
compatible with experimental data. 

As an additional information, MiniBooNE also gives a ratio more directly related to 
the measurements, namely the ratio of pion-like events (defined as events with exactly one 
\T and one tt + escaping the struck nucleus) and quasi-elastic signal (defined as those with 
one \i~ and no pions). In our language the last quantity represents the total Np — Nh 
(N = 1, 2, 3, including the quasi-elastic for = 1) exclusive channel. We have compared 
this second experimental information to the ratio between our calculated pion production 
(which however ignores final state interactions) and our total Np — Nh contribution to 
the total charged current neutrino cross-section (lower panel of Fig JT9|) . There is an 
appreciable difference between the two curves of Fig]T9] the one in the lower panel is 
reduced due to a large 2p—2h component in the Np—Nh cross-section, which increases the 
denominator. The comparison with the experiment shows an agreement up to E v ~ 1.2 
GeV. Final state interactions for the pion, which are not included, are expected to reduce 
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our result at the level of 15 %, still maintaining an agreement. 

A new result has been presented at Nulnt09 by SciBooNE [16|. It is the ratio of the 
total neutral current 7r° production cross-section to the total charged current cross-section 
at the mean neutrino neutrino energy of 1.16 GeV. They obtain the preliminary value: 

(7.7 ± 0.5(stat.)lg;t(sys.)) ■ 10~ 2 . (14) 



a(NCn ) n , Kr;ir )+(Uf i; 1A _ 2 

(T[L>L" 



TOT 



Our prediction for this quantity, including coherent contribution and a factor 2/3 for NC 
incoherent pion production to single out n° contribution is: 

a(NCn ) _ in _ 2 



a(CC- 



TOT) 



7.9 -10 , (15) 



which fully agrees with data. 

A general comment on the comparison with data: nearly all the ratios that have 
been discussed, except the final-state-interaction-corrected MiniBooNE result of the upper 
panel of FigHH are sensitive to the presence of the np — nh (n = 2, 3) component in the 
cross-section. Since the size magnitude is not so well tested, we can investigate what 
becomes the comparison with data in the extreme situation when we totally suppress this 
contribution. For the last ratio discussed we obtain 

a{NC no) 9.8 ■ ID" 2 , (16) 



((j(CC TOT ) - o{CC np _ nh )) 



appreciably above the experimental value. 

As for the SciBooNE upper limit of the ratio of the 7r + coherent to total charged 
current cross-section, our prediction at E u =l.l GeV, which was 0.71-10 -2 , without np—nh 
becomes 0.89-10 -2 , further above the experimental bound of 0.67-10 -2 . 



4.3 Quasi-elastic cross-section 

A new preliminary result on absolute cross-sections has been presented by the MiniBooNE 
collaboration This group gives in particular the absolute value of the cross-section 
for "quasi-elastic" events, averaged over the neutrino flux and as a function of neutrino 
energy. The comparison of these results with a prediction based on the relativistic Fermi 
gas model using the standard value of the axial cut-off mass Ma = 1.03 GeV/c 2 reveals a 
substantial discrepancy. In the same model a modification of the axial cut-off mass from 
the standard value to the larger value Ma = 1.35 GeV/c 2 is needed to account for data. 
A similar conclusion holds for the Q 2 distribution The introduction of a realistic 



spectral function for the nucleon does not alter this conclusion [32| . 

As a possible interpretation we question here the real definition of quasi-elastic events. 
As already discussed above, the nuclear medium is not a gas of independent nucleons, 
correlated only by the Pauli principle, but there are additional correlations. The ejection 
of a single nucleon (denoted as a genuine quasi-elastic event) is only one possibility, and 
one must in addition consider events involving a correlated nucleon pair from which the 
partner nucleon is also ejected. This leads to the excitation of 2 particle-2 hole (2p — 2h) 
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Figure 19: Ratio of the z/^-induced charged current one n + production to quasi-elastic 
cross-section as a function of neutrino energy. 
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Figure 20: "Quasi-elastic" u^- C cross-section per neutron as a function of neutrino en- 
ergy. Dashed curve: pure quasi-elastic (lp — Ih) cross-section; solid curve: with the 
inclusion of np — nh component. The experimental MiniBooNE points are taken from 
H. 



states which have been abundantly discussed throughout this work. In the spin-isospin 
channel the correlations, mostly the tensor ones, add Ip — lh strength to the lp — lh events 



661 ]. At present, in neutrino reactions, such events cannot be experimentally distinguished 
from the genuine quasi-elastic events and must be considered simultaneously. Notice that 
the standard lower value of the axial mass, M4=1.03 GeV/c 2 , results from deuterium 
bubble chamber experiments. In this case the effect of tensor correlation is also present 
but at a lower level since deuteron is a dilute system. Our sum of the combined 12 C 
quasi-elastic cross-section and the 2p — 2h one is displayed in Figj20l This prediction fits 
the experimental data excellently, better than expected in view of the uncertainties of our 
lp — lh cross-section. As for the flux averaged "quasi-elastic" cross-section per neutron the 
experimental value is 9.4 10 -39 cm 2 (with a normalization error of 11%). Our prediction 
for this quantity is 6.3 10 -39 cm 2 without 2p — 2h contribution and 9.0 10~ 39 cm 2 including 
it, a value more in touch with the experimental one. 

In view of the importance of the issue we want to investigate if this large 2p — 2h 
contribution is a genuine effect and not an artifact of the particular parametrization that 
we have used for the bare 2p — 2h channel. For this, we introduce a different approach 
which exploits a microscopic evaluation by Alberico et al. 66] of the 2p — 2h contribution 
to the transverse magnetic response of (e, e') scattering. It does not have the shortcomings 
of our previous parameterizations which have no momentum dependence. In the previous 
case the maximum of the 2p — 2h response R^^h always lies at low energy, u ~ 50 MeV, 
irrespective of the momentum, separating at large momentum from the quasi-elastic peak 
which instead gets shifted at larger energies. A similar feature exists in the NA part. 
This is not realistic and below we sketch a possible way for improvements. The aim is 



to extract the 2p — 2h responses from the results of Alberico et al. [6_6j, although they 



23 




CO [GeV] 

Figure 21: Comparison between the CC - 12 C np — nh (n — 2,3) differential cross- 
sections deduced from the two different parameterizations of the 2p — 2h components. 
Thin lines: parametrization of Sec J3.3l and used throughout the whole article (denoted 
"old"); bold lines: parametrization of Sec J4.3l ("new"). 



are available for a limited set of momenta and energies and they concern iron instead of 
carbon. We have thus performed extrapolations both to cover all the kinematical region 
of neutrino reactions and to go to the 12 C case. For the set of R aT tT\{uj , q) values that 



we could extract [66j we have observed an approximate scaling behavior with respect 

2 2 

to the variable x = ■ A parametrization of the responses in terms of this variable 
allows the extrapolation needed to cover the full neutrino kinematical region and we have 
now the new responses, R^p-^hi 10 S Q) an d -^2^-2/1(^5 <?) m a ^ the range. For the AA part, 
which is not well covered in [66] we have kept the previous parametrization, which already 
presents a proper q dependence owing to the contribution of the in-medium A width [70|- 



Another remark is in order. The evaluation of Ref. 66j of the 2p — 2h channel does not 



reproduce pion absorption in nuclei at threshold, as observed by the authors. It gives a too 



large value for the absorptive p-wave optical potential parameter [75|], ImCo — 0.18m~ , 
instead of the best fit value ImC — 0.11m" 6 . To be as consistent as possible with our 
previous parametrization, which comes from pion absorption, we have applied to our 
scaling function the reduction factor j^§. The nuclear mass dependence is taken care of 
with the introduction of the Levinger factor, L, which fixes the number of quasi-deuteron 
pairs in the nucleus defined as L ZN/A. We rescale the iron results by a factor r, ratio 
of the Levinger factors, for the two nuclei. It is r = 0.8 according to the A dependence 
of the Laget formula [83j or a similar value, r ~ 0.75 from [84j. Altogether the global 
reduction factor applied to the iron scaling function is ~ 0.5. 

Since in the previous case the RPA have little effect on the 2p — 2h component, we 
introduce directly the bare new 2p — 2h quantities in the neutrino cross-section. The 
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Figure 22: Same as FigEZUl including also our curve (QE + np — nh "new") with the new 
parametrization for the 2p — 2h component (bold line). 



influence of the new modelization of the 2p — 2h is displayed in Fig. [21] where the bare 
partial and total np — nh differential CC neutrino cross-sections at E v = 0.7 GeV are 
shown both for the previous parametrization and for the new one. The energy behaviors 
are quite different, the NN contribution is no longer localized at small energies but is 
spread over the whole energy range, a more realistic feature. A similar behavior occurs for 
iVA part. However the integral over the energy, a np - n h(E u ) , is practically not modified. 
As a consequence adding this contributions to the quasi-elastic cross-section we reach a 
similar conclusion as before about the important role of the multi-nucleon channel, as 
illustrated in Fig. El 

It indicates that, in the nuclear medium, neutrino in this energy range do not interact 
only with individual nucleons but also with pairs of nucleons, mostly n—p pairs correlated 
by the tensor interaction. The spin dependent part of the neutrino interaction with such 
a pair is stronger than with the same two nucleons when isolated. This increase manifests 
itself trough the 2p — 2h strength which adds to the lp — lh part, an effect simulated 
by an increase of the axial cut-off mass. Quantitatively a confirmation on the theoretical 
side of the exact magnitude through a detailed microscopic calculation of the bare 2p — 2h 
response, which will then be inserted in our RPA formalism, would be helpful. Also an 
experimental identification of the final state would be of a great importance to clarify 
this point. In particular the charge of the ejected nucleons will be quite significant. Since 
tensor correlations involve n—p pairs, the ejected pair is predominantly p — p for charged 
current and n—p for neutral current. This predominance has the same origin as for p-wave 



7r absorption by nuclei where n — n emission is favored over n — p emission [72j, [79 



5 Summary and Conclusion 

We have studied neutrino interactions with light nuclei which enter the targets of present 
or future experiments. Our theoretical tool is the theory of the nuclear response treated 
in the random phase approximation (RPA), a well established technique for the treat- 
ment of electromagnetic or weak interactions with nuclei and which have been used also 
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for strongly interacting probes. The crucial element of the RPA treatment is the p — h 
interaction, in particular for the spin-isospin one, which has been taken from the accu- 
mulated knowledge on the spin-isospin responses. The main merit of this approach is to 
allow unified description of various channels. It has some limitations which restrict the 
energy range of the neutrino to a region below ~ 1.2 GeV. For instance the only nucleonic 
resonance incorporated in the description is the Delta resonance. Multi-pion production is 
also ignored, as well as most meson exchange effects. Moreover, although both the Delta 
propagator and the center-of-mass correction are the relativistic one, not all relativistic 
effects are included in a systematic way. 

The final states considered are the quasi-elastic, 2p — 2h, 3p — 3h ones, and coherent or 
incoherent pion production. Some channels have the problem that final state interactions 
are not incorporated. This is the case for incoherent pion emission where the produced 
pion can be absorbed on its way out of the nucleus leading to a multi-nucleon state. 
Incoherent pion production is therefore overestimated and the multi-nucleon channel ac- 
cordingly underestimated. This effect is visible in the scattering of physical pions on 12 C 
in the region of the Delta peak. For a light nucleus such as 12 C the effect is limited but 
it becomes more serious in heavier nuclei. Our method should be supplemented by an 



evaluation of the final state interaction, for instance by a Monte Carlo method 61 
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The coherent channel is particularly interesting although it represents only a small 
fraction of the total pion emission. It does not suffer from the previous limitations as 
final state interactions are automatically incorporated in the RPA treatment which is 
particularly suited for this channel. Moreover, it is the only channel which is dominated 
by the isospin spin-longitudinal response where collective effects are very pronounced 
while they remain moderate in the other channels. The difference between the first order 
term with one bubble (with Delta excitation) and the full RPA chain is quite large. In 
this context we have used as a test of our spin longitudinal response the experimental 
data on elastic pion scattering in the Delta region. It offers a direct test of the forward 
coherent neutrino cross-section to which it is linked through the Adler theorem. Except 
for low pion energies near threshold (u < 200 MeV) where Adler's theorem fails, the 
elastic cross-section can be used to extract the forward neutrino coherent cross-section as 
in the method of Rein-Sehgal. 

For the evolution of the partial cross-sections with mass number in order to reach 
the 40 Ca region, our description indicates that, apart from the coherent pion production 
which evolves differently, the other partial cross-sections scale essentially as the nucleon 
number. Final state interactions obviously will destroy this scaling. 

We have compared our predictions with the available experimental data. Our ratio 
for the coherent pion production over the total neutrino cross-section is just compatible 
with the experimental upper limit. Another test concerns the ratio for charged currents 
7r + production to the quasi-elastic cross-section. A delicate point in the experiments 
concerns the definition of a quasi-elastic process and its separation from np — nh which 
the experiment does not distinguish. In one set of data a correction is applied to obtain 
a genuine quasi-elastic cross-section and it is corrected as well for final state interaction. 
In another set of results a generalized quasi-elastic is introduced, defined as events with 
only one lepton. In this case our 2p — 2h and 3p — 3h should be added to the quasi- 
elastic component. Both lead to successful comparisons with the two sets of experimental 
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data. Further data involve a ratio of neutral current n production to the total neutrino 
cross-section for charged currents. Here again our evaluation agrees with data. It is quite 
encouraging that the comparison with present experimental data is essentially successful. 

A distinct feature of our approach, and one of our significant results, is the large 2p—2h 
component. It affects all the measured ratios discussed in this work. At the present level 
of accuracy we have not found in these ratios any contradiction to its presence. It is 
also supported by preliminary data on the absolute neutrino quasi-elastic cross-section on 
carbon. We suggest that the proposed increase of the axial mass from the standard value 
to a larger one to account for the quasi-elastic data, reflects the presence of a polarization 
cloud, mostly due to tensor interaction, which surrounds a nucleon in the nuclear medium. 
It translates into a final state with ejection of two nucleons, which in the present stage 
of the experiments is indistinguishable from the quasi elastic final state. Although the 
existence of such 2p — 2h component is not in question, for a fully quantitative evaluation 
we plan to improve the description of the multi-nucleon final states by a microscopic 
treatment. Future precision experiments, such as T2K, will be able to identify final 
states, namely p — p pairs for charged current and n — p pairs for neutral current, and 
bring an experimental elucidation of this intriguing effect. 
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A Inclusive neutrino-nucleus cross-section 

The invariant amplitude for the lepton-nucleus cross-section, Eq. (DO), results from the 
contraction between the leptonic L and the hadronic H tensors 

\T\ 2 = L 00 W 00 + L 33 W 33 + (L 03 + L 30 )W 03 + 

(Lu + L 22 )W U ± (L 12 - L 2l )W 12 \ + ^\ . (17) 

The various L are the component of the leptonic tensor 

Lfj,u = 8(k^k v + k v k ^ g^k.k =p i^^vapk k ) (18) 
while the W of the hadronic one 




where E q N = (q 2 + M N 2 f/ 2 and E\ = (q 2 + M A 2 ) 1/2 . This decomposition takes into 
account the different channels of particle-hole excitations. 
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The various leptonic tensor components are: 



L o = S(kok'o + kk' cos 9), 

L 33 = 8(2k 3 k' 3 + k k' - kk! cos 0), 

L 03 + L 30 = 16(k k' 3 + k 3 k' ), 

L n + L 22 = 16(2kl + k k' -kk' cos 0), 

L12 - L21 = -i 8(/c A;'3 - k 3 k' ), 

with k 3 = —(k' cos9 — k), 

k' 3 = —{k' - kcos0), 
Q 

, , , sin # 
and ki = k 2 = 20 

For the hadronic tensor components we keep only the leading terms in the development 
of the hadronic current in p/M, where p denotes the initial nucleon momentum. Marteau 
investigated the importance of the momentum terms and found them to be small. The 
components are related to the various nuclear responses R as follows 
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(21) 



For sake of illustration we give the explicit expression of H 00 : 



rOO 
Ipp, 



H 00 = J2 H " 

P,P'=N,A 

t?NN , nO nO r>NN , q ,q0 ,o0 r>NA , nO aO nAA 



a N a N R T + P N P N R l + 2(3 N (3 A R l + j3 A f3 A R l 



(22) 
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The quantities a, (3, 7 and 5 are expressed in terms of the usual form factors, namely 
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(23) 



We have introduced in the time component of the axial current a renormalization factor 
G* A = Ga(1 + 5) to account meson exchange effects which are known to be important in 
this channel 68]. Even with the large value 5 = 0.5 the effect of this renormalization is 



small. The most affected channel is the coherent, which is reduced by ~ 10%. 

A.l Spin longitudinal contribution to the quasi elastic cross-section 

We consider the limit of vanishing lepton mass. In this case the relevant leptonic tensor 
components reduce to 



^00 



4[(* + k'Y - q'\ = \L 



OJ' 



33 



1 q 

"o — (A)3 + -^30 )• 

2 uj 



(24) 



The longitudinal contribution to the quantity |T| 2 involves 



0n -^00 + Pn L 33 + P%Pff(L Q0 + L33) — NIj 2 G 2 a Lqq 



(E q N + M N y + 7 



qE q N + M N 



(25) 

Neglecting the struck nucleon momentum, the transferred energy oj in a quasi-elastic 
process is oj = E q N — M N , which implies the bracket on the r.h.s. of Eq. (T25l) to vanish. 



B Particle-hole polarization propagators 



B.l Bare 



In this Appendix we give the expressions of the bare particle-hole polarization propagators. 
The nucleon-hole polarization propagator is the standard Lindhard function [86J. 
For the Delta-hole polarization propagator we use the relativistic expression 



n A -h(g) 



32M A f d 3 k 



9 



(27T) 



; 6{k F - k) 
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s-Ml + iM A T A 
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(26) 
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where s and u are the Mandelstam variables. Ma = Ma + 4Q(MeV)j^ is the mass of the 



A in the nuclear medium and Ta is the in medium Delta width. The last two quantities 
are taken from 



70 



For the 2p — 2h polarization propagators we consider only the imaginary parts. Their 
expressions, which represent an extrapolation of threshold results of [72| are 



Im(U° NN ) 
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(2M N + m n ) 2 
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C 3 $ 3 (a;) 



(uj + M A - Mv) 2 



(27) 



The Cj constants are set to C x = 0.045, C 2 = 0.08, C 3 = 0.06, while the $j(o;) include 
phase space, pion and rho propagators. 

B.2 RPA 

Here we define the RPA expressions of the response functions for finite nuclei. 

First we introduce the projection of the bare propagators on the Legendre's polyno- 
mials Pl through 



U^ L \uj,q,q') = 2tt / duP L (u)U°(uj,q,q'), 



k ' F ( R) {uJ,q,q 



2tt / duP L (u)U° k 



(28) 



where q = \q\, q' = \q'\, u = cos(q,q'). 

Starting from Eqs.([5]) and f j28~l) . after some algebraic manipulations, one obtains 



n°( L \u,q,q') = 47T^(2ii + l)(2l 2 + 1) ( 

x J dRR 2 3 h (qR) Jh (q'R)T]: 



h h L 




0(h) (t, a n> 



(29) 



with the usual three- j symbol and /-order Bessel function ji(x). 

This is the starting point for the calculations of isovector and spin-isospin response func- 
tions. 

The free isovector (or charge) response function can be expressed through 
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Vv^2J+l 
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' 7T *yi 4vr 
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(30) 
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The RPA isovector response function 

RZ N (^Q) = -^Im[n ccjvJV (^,g,q),] = -^^17^ Im [ n S>>?>?)] , (31) 
is obtained solving the following equation 

n£>> ?, g') = n»(", g, g') + / ^ nj^o,, g , fc) ^(*) n^( W) fc, </)■ (32) 

For the spin-isospin longitudinal and transverse responses we introduce the following 
quantities 
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27TT forL=J-l, 



>.JL 



J 



2J+1 

1 for L=J 



for L=J+1, (34) 



Note that, in general, for finite systems 7^ 0. 

The bare responses in a particular channel k (k — QE, 2p — 2h, ...) are given by 

j 

with x,y = l,t, referred to the longitudinal or transverse channel, and PP' = N, A. 
The second term of Eq.(TT0l) in the channel k, namely 

Imn (fc) = |1 + IIVf Imn^, (36) 
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with II the full polarization propagator, explicitly writes 



+ 



(2vr) 3 

^ ' QR 



WW 1 



dpp 2 dp' p' 2 



(2tt) 3 (2tt) 3 



'HH' 



n «i^, ^ p') (pO n -'^P' p'> <?') > (37) 

where x, y, w;, w', z, z' = I or t and P, P', Q, Q', R,R' = N,A. 

The solution of this equation leads to the corresponding response functions 

9) = ~ E ^ Im [ n SL, (<". «)]■ ( 38 ) 
j 

In our calculations the maximum multipole number is set to J = 25 which turns out to 
be sufficient to reach the convergence. 

The first term of Eq.lTTOl). which represents coherent processes, explicitly writes 

n (i^>>^') = / ^(n&(^?,p))^(^(p))nJ2?p,( WlPig 



( n Lp7v( w ; 9> ns; P ,(o;, g^, g') 



167T 2 ^ PiVV ^>H,Hn)J ^lyNP' 

+ ^4" ( n 2w( w > 9, 9tt))* n[^p,(w, g WJ g 7 ) 
+ ~7n2 ( n ifpA(^ ; 9, 9tt))* n[^ P ,(w, q^, q') 

c*2 



+ £f (nS A (^, 9, nf v lp,(u;, q ni q')) , (39) 



where 



Q 2 



Im (T4) = Im C n - ^ 2 ■ , = -i C w 7T qr^(g 2 - <£) = -i tt ^ <J(|g| - 

\ uj z — q z — m£ + zr] J 2 

(40) 

with the generic Nucleon- or Delta-pion coupling constant and q^ = y 1 'uj 2 — m 2 . 
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